Review of Applied Physics Vol. 2 Iss. 2, June 2013 


ww w. seipub . or g / rap 


Analysis of the Behaviour of an Optimized 
Magnetorheological Brake 

Giuseppe Marannano 1 ; Gabriele Virzi' Mariotti 2 ; Cedomir Duboka 3 

12 Dipartimento di Ingegneria Chimica, Gestionale, Informatica, Meccanica Universita di Palermo, Viale delle 
Scienze 90128 Palermo, Italy, 3 University of Belgrade, Faculty of Mechanical Engineering, SRB-11120 Belgrade 35, 
K. Marije 16, Serbia 

n giuseppe. marannano@unipa.it; 2 gabriele.virzimariotti@unipa.it; 3 cduboka@eunet.rs 


Abstract 

This paper presents the optimization process of MR - 
magnetorheological brake in order to define its configuration, 
and to reach the requested vehicle braking torque. FEM 
analysis previously carried out, concerning structural sizing 
of such a brake, showed that both braking torque and brake 
mass did not satisfy the requirements. Due to the necessity 
of limitation on the suspended mass of a motor vehicle, the 
optimization of the form of MR brake stator has been 
executed, including determination of its contribution to the 
total mass of the brake. Thus a percentage reduction of 
approximately 40% was obtained. In order to estimate the 
brake temperature which results from the transformation 
into heat of a vehicle kinetic energy during brake application, 
a thermal Finite Element Analysis is carried out. In 
particular a fade braking test composed of snub braking, i.e. 
repeated cycles of acceleration and consecutive braking 
applications is performed. 

Keywords 

Magnetorheological Brake; FEM Optimization; Magnetic FEM; 
Thermal FEM; Fade Braking Test 

Introduction 

Sizing of MR brake was executed in the preliminary 
work (Marannano, Virzi Mariotti and Duboka, 2011), 
following the criteria available in the literature (York, 
1997; Carlson LeRoy, Holzheimer, Prindle, Marjoram., 
1998). The requirements (Carlson, 1999; Brauer, 2006) 
for an automotive brake were not found to be satisfied 
in this way. In fact FEM Analysis, performed on the 
model of MR brake, reveals the possibility to bring 
modifications to the magnetic circuit, with the aim to 
improve the requested performances, in terms of the 
maximum braking torque value and the minimal 
weight of the device. 

The procedure of MR brake design is studied, and 
several papers are found in the literature (Karakoc, 
Park and Suleman, 2008; Farjoud, Vahdati and Fah, 
2008; Yang, Huang and Kang, 2007; Li and Du, 2003; 


Liu, Li, Kosasih and Zhang, 2006; Park, Stoikov, da 
Luz and Suleman, 2006; Sukhwani and Hirani, 2008; 
Zhou, Chew and Hong, 2007). Several modifications 
on the original procedure have been proposed, but the 
FEM Analysis showed that the analytical procedure 
did not enable the values of the required braking 
torque to be obtained. Additional optimization is 
carried out, to obtain the achievement of requirements 
(Yang, Duan and Eriksson, 2008; Jaindl, Kostinger, 
Magele and Renhar, 2007; Park. Luz and Suleman, 
2008). 

Compared to similar analyses found in literature, MR 
brake in present research consists of two discs instead 
of one (Nguyen and Choi, 2011); in order to contribute 
to the feasibility of the brake design. Magnetic circuit 
optimization process is executed in two phases, using 
parametric FEM model of the brake. The primary 
objective of the first phase is the determination of the 
optimal geometric parameters of the magnetic circuit, 
to obtain values of braking torque not lower than the 
nominal (or designed) value - Mfnom, representing 
the limit of a typical motor vehicle brake. Starting 
from the obtained configuration, the second phase of 
optimization was deemed to search an optimal form of 
the stator in ferromagnetic steel that is the novelty 
regard to previous works. The decreasing trend of the 
magnetic induction B versus the decreasing distance 
from the rotation axis suggests a shape of variable 
thickness of the stator that is modeled by a spline, 
while the coordinates of the control points are chosen 
as design variables; thickness of MR fluid (4 = 1 mm) is 
not modified during the optimization procedures. 
Thus a sensitive reduction of the mass of the entire MR 
brake may be obtained, without significant reduction 
of the braking torque performance, to acquire very 
original shape of the brake geometry, with 
simultaneous achievement of requirements in terms of 
both brake torque and mass. 
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Finally the fade behaviour of the brake is analysed 
from the viewpoint of dynamic thermal analysis, 
to simulate a braking test with the maximum 
temperatures that significantly affect the efficacy 
of MR fluid (Nguyen and Choi, 2011, Park, Luz 
and Suleman, 2008; Wiehe, Noack and Maas, 
2007). 


Optimization process 

The process of optimization (Kodiyalam and Saxena, 
1994) of an element consists of search on a whole of 
design variables (state variables), minimum value 
assumed by a determined objective function, and 
satisfication of the requirements demanded by the 
system. The mathematical formulation of a typical 
process of optimization is: 

Minimizing f ({x}) 

Subjected to: gj ({x}) <0 j = 1,2, . . . m (1) 


being: 


X i , min ^ X, < X- max 

{x} = [x,,x 2 ...xj 


i = 1,2, /i 


where /*({x}) represents the objective function, {x} is 
the vector of independent variables of design, in 
number equal to n, subjected to constraints imposed 
by k.m.n’-Vmax J while &({*}) represents the whole 

of state variables, in equal number to ;, functions of the 
design variables that the system must satisfy. 

Therefore the optimization represents a procedure 
pointing out the best design variables, inside an 
admissible domain, that render minimum (or 
maximum) a determined function of the variables 
governing the problem. 

The search for "an optimal" solution to the problem 
can be take advantage of a great variety of methods. 
The choice of the method is a direct consequence of the 
type of problem under investigation. 

A distinction is done between problems of conditioned 
type, in the admissible design configuration with 
respect of imposed conditions (constraints), and 
problems of unconditioned type that have not 
restrictions. This last have simpler solution but rare 
presence in engineering (Nocedal and Wright, 2006), 
(Elster, 1993). However, several methods may be used 
to transform a problem of conditioned type in 
unconditioned one (method of the penalty functions). 

Another important factor concerns the existence of 
eventual local minimums of the objective function. 


Such type of problem can be resolved by means of 
various methods, among which the execution of 
several analyses are started by a whole of design 
variables inside the permissible domain. A complete 
description of the state variables is presented in the 
next paragraphs 3.1 and 3.2. 

1) Ansys optimization 

ANSYS executes a series of cycles of analysis- 
appraisal-modification of the model testing in an 
iterative way the design variable (DV), state 
variable (SV) and the objective function (OBJ). Such 
variables are defined in parametric way, by means 
of language AD PL. The constraints on both the 
design and the state variables are imposed by the 
user. Their choice, with that of the objective 
function, represents one of the most important 
elements of the decision phases to really obtain a 
solution that represents the optimal design. A 
whole of values of design variables satisfying the 
imposed constraints for the problem under 
investigation, is indicated like feasible by the 
software; otherwise it is indicated as infeasible and 
eliminated by the whole of the solutions. The 
optimization procedure uses the method named 
Subproblem approximation method, for the 
minimization of the objective function. The 
software uses some techniques of optimization, by 
means of which the efficiency of the process is 
increased and the design random generation is used. 
The better solutions constituting the point of 
departure for the successive procedure of 
optimization, may be selected by means of the 
instruction OPSEL. 

First phase of optimization of the magnetic circuit 

The first step of the optimization process of the 
magnetic circuit regards the choice of the objective 
function and the design variables, with the 
determination of the constraints (variables of state). 

The primary objective is the determination of the best 
geometry of the magnetic circuit to obtain a value of 
braking torque equal to requested nominal value Mfnom. 
An objective function must be imposed that considers 
feasible the geometrical configurations of the circuit, for 
which the condition M f > M j- nom is fulfilled, being Mf 

the actual braking torque. 

The objective function is defined in the form: 

f( x ) = K- M f (2) 
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The term K in (2) represents a constant, introduced 
with the aim to avoid numerical problems, due to 
negatives values of fix) during the optimization 
process. A value much greater than the braking torque, 
equal to 10 4 , is chosen. 

An important aspect is the mass of the brake, which 
has to assume to be the lowest possible values, 
because too heavy brake greatly increases the overall 
weight of the vehicle. Since the total mass of the brake 
is mainly determined by the size of the magnetic 
circuit (greater density of ferromagnetic steel and of 
the copper of the coil), a maximum acceptable limit 
has to be considered, by introducing the state variable: 


Value of Pmax fails to reduce the search domain of the 
optimal solution. Referring to the mass value of the 
model, obtained from a first approximate sizing, a 
maximum value of Pmax=0,8 P= 33,8 kg is imposed. 
The ratio 0,8 is not exclusive and the objective is a 
lower mass than the initial one. If a lower value of the 
ratio is chosen, the same value of the mass is obtained 
as the final result. 

Mass of the MR brake can not be compared to the 
hydraulically or mechanically operated braking 
systems which are characterized with the presence of 
additional components (brakes pump, oil tank, 
hydraulic circuit, mechanical linkage, etc.) replaced by 
electrical components in the case of MR brake. 
Another constraint is the upper limit of the magnetic 
field B in the stator steel. Indeed, if the value of 
magnetic field in steel CIO is close to the saturation 
value, the circuit is unable to produce an increase in 
the magnetization of the MR fluid, increasing the 
current supplied to the system. This causes loss of 
linearity in the electromechanical characteristic 
(torque-current) of the device. A state variable is then 
defined taking into account the saturation of steel, 
assuming the following expression: 

B<0,9B sat (4) 


TABLE 1 DESIGN VARIABLES AND RELATIVE CONSTRAINTS. 


Design variables 

Range of variation [mm] 

Inner radius of active surface n 

40-140 

Radial thickness of the coil Lc 

5-20 

Disk thickness Sd 

2-10 

Stator radial thickness SI 

5-30 


The choice of the limits within which the design 
variables change takes into account two factors. On 
one hand too large range do not adequately represent 
the entire admissible domain, while on the other hand 


too narrow range may exclude the best solutions. The 
selected design variables and their range of variation 
are listed in Table 1. 

TABLE 2 RESULTS OF THE FIRST PHASE OF OPTIMIZATION. 


Objective function 

Obtained values 

Braking torque [Nm] 

1055,34 

Mass [kg] 

29,75 

Design Variables 

Optimal values [mm] 

Inner radius of active surface n 

133,5 

Coil radial thickness Lc 

9,2 

Disk thickness Sd 

3,2 

Stator radial thickness Si 

15,34 



FIGURE 1 MAGNETIC FIELD DISTRIBUTION. 



1) Results 

Table 2 shows the results of the first phase of 
optimization of the magnetic circuit, where the best 
configuration is reported. Analysis of Table 2 
shows that, in the best configuration, the MR brake 
is capable to generate a value of braking torque 
above the minimum required. The mass of the 
brake is lower than the preset limit and the 
reduction is approximately 30% compared with the 
value obtained from the preliminary analysis. 
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Figure 1 shows the distribution of magnetic field in 
the circuit in the optimal solution. The maximum 
value of Bf at the fluid gap is now equal 
to B f = 0,8 T , which is sufficient to obtain the 

required magnetization. Fig. 2 shows the flow lines 
of magnetic field. 



FIGURE 3 MAGNETIC FIELD INTENSITY H. 
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FIGURE 4 BRAKING TORQUE VERSUS NI FOR EACH 
ACTIVE SURFACE (FIRST OPTIMIZATION). 


Figure 3 shows the map of the intensity of magnetic 
field H in the circuit. The mean value of magnetic 
intensity is equal to Hf = 140 kA/m, corresponding 
to a shear stress tH = 52 kPa, which is very close to 
the maximum shear stress of MR fluid. Figure 4 
shows the load variation versus the magnetomotive 
force (NI) supplied to the system, obtained from 
interpolation of data for each active surface. An 
almost identical trend can be observed for each pair 
of external (gap 1 and 4) and internal (gaps 2 and 3) 
gaps. The contribution of each gap to the value of 
total braking torque is slightly higher for each gap 
on the outside, near to the maximum value of 
magnetomotive force. The electromechanical 
characteristic of the brake may be estimated adding 
the contributions of braking torque for each gap 
fluid, as Figure 5 shows. 



NI [A] 

FIGURE 5 ELECTROMECHANICAL CHARACTERISTIC OF 
THE BRAKE (FIRST OPTIMIZATION). 

The torque has a roughly linear trend to vary NI, 
except the last section where the curve decreases its 
slope. This completely predictable behaviour is due 
to values of current near to the maximum value, so 
that the magnetic field in the circuit tends to reach 
the saturation value of ferromagnetic steel. 

Braking torque assumes the value to be zero in the 
absence of applied current; due to the fact that the 
contribution of the Newtonian term of the shear 
stress is neglected in the calculation of total braking 
torque. 

Torque Mf n due to viscosity in lack of magnetic 
field is given by: 

( 5 ) 

Where rj = 0,280 Pa is the viscosity, E, = 0,001 m is 
the gap, co is the angular velocity (which is equal to 
118,2 rad/s for initial speed of 120 km/h and the 
tyre rolling radius of 0,282 m), re and n are external 
and internal radius of the disk respectively. A 
value of the torque equal to 33 Nm is thus obtained; 
that is 3,1% of the value reported in table 2. 

Second phase of magnetic circuit optimization 

Figure 1 shows that the stator in ferromagnetic steel is 
characterized with regions at low density of magnetic 
flux. In particular, the value of B decreases in radial 
direction towards the regions near to the axis of 
rotation of the brake. For the principle of continuity of 
magnetic flux, considering two generic sections of the 
stator, one has: 

fa=t 2 =B l -A 1 =B 2 -A 2 ( 6 ) 
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Since the value of magnetic field B is inversely 
proportional to the size of the section crossed by the 
flow, the thickness of the stator can be reduced in axial 
direction in zones where B is low, compared to the 
magnetic saturation of the steel. 

In light of these considerations, the second phase of 
the optimization process is useful to find a suitable 
form of the stator, to reduce the mass without affecting 
the system performance significantly. 

1) Changes to FEM model 

A modification of the geometry of the model is 
necessary to search an optimal form of the brake 
stator. The stator is achieved in a completely 
parameterized way by means of bottom up 
procedure. The outer contour is achieved by the 
use of Spline Curves, passing through checkpoints 
set. This approach allows obtaining an external 
profile of the stator, through the variation of the 
vertical coordinate of these points. The procedure 
has already been described in the previous 
paragraph, so that the intermediate steps are the 
same; the optimization parameters are indicated in 
the next paragraph. 

The active surfaces of the model are implemented 
in parametric form by direct generation of nodes 
and elements for the subsequent calculation of the 
braking torque. Figure 6 shows the geometry of the 
brake, the control points and the stator thickness. 

2) Definition of the optimization parameters 

Since the second phase of optimization is devoted 
to the research of the lower mass of the brake, the 
objective function is represented by this variable: 

m=p < 7 ) 


Changing the shape of the stator as well leads to an 
increase in the value of the magnetic field passing 
through it. This leaves an increase of the total 
reluctance of the circuit, reducing the intensity of 
magnetization in MR fluid. The introduction of the 
new state variable is necessary. 


M f >M 


fnom 


( 8 ) 


to avoid the reduction of the performance 
requirements in terms of torque produced by the 
brake. The number of design variables is greater 
compared to the previous optimization. The 
establishment of the stator thickness at control 
points requires some essential considerations, 
because their variation during the optimization 


process cannot be arbitrary, to avoid geometric 
configurations that may be difficult to implement. 
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FIGURE 6 CONTROL POINTS IN THE BRAKE GEOMETRY. 
TABLE 3 DESIGN VARIABLES AND CONSTRAINTS. 



FIGURE 7 DISTRIBUTION OF THE MAGNETIC FIELD B IN 
THE SECOND OPTIMAL CONFIGURATION. 

The decreasing trend of the magnetic field through 
the stator (Fig. 1) suggests, like possible 
configuration, a stator thickness gradually growing 
in radial direction away from the axis of rotation. 
On that basis, a dependency relationship between 
the control thicknesses from ti to U may be 
established, in order to obtain really feasible results 
of the optimization. In particular, for given value of 
ti, the following relationship is imposed: 

t M =t,+ A, / = !... 3 (9) 
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Ai represents the successive increase in thickness, 
giving a tapered form towards the outside of the 
stator. Table 3 shows the selected design variables 
and their ranges of variation established according 
to the results obtained from the first phase of 
optimization. 

Maximum radius of the brake is Rmax = 200 mm 
and thickness of MR fluid is ^ = 1 mm; excluding 
from modification during optimization procedure. 

3) Results 

Table 4 shows the results obtained for the geometry 
optimization. The second phase of optimization 
allows a reduction of brake mass of about 14% 
compared to the value obtained in the first phase. It 
is also important to note that, despite the increase 
of the reluctance of magnetic circuit, due to the 
increase of the magnetic field B, the value of 
braking torque is reduced by less than 3,5%. Figure 
7 shows the distribution of the magnetic field in the 
circuit, in the plane. The maximum value of B in 
the stator is greater than the previous, but lower 
than the limit of saturation of the steel. 


TABLE 4 RESULTS OF THE SECOND PHASE OF OPTIMIZATION. 


Braking torque [Nm] 

1020 

Mass [kg] 

25,8 

Design variables 

Optimum values [mm] 

Inner radius of active surface n 

126,8 

Coil radial thickness Lb 

10 

Disk thickness Sd 

2,5 

Stator radial thickness Si 

14,6 

Stator axial thickness ti 

12,2 

Stator axial thickness ti 

12,8 

Stator axial thickness h 

17,1 

Stator axial thickness U 

17,1 

Stator axial thickness ts 

16,2 

Fillet radius r s 

10 


Figures 8 and 9 show the electromechanical 
characteristic of the new configuration, both for 
individual disks and the entire brake, in a similar 
way to the previous optimization phase. 

Thermal analysis of the MR Brake 

Working principle of the MR brake foresees the 
presence of sliding between the active surfaces in 
relative motion. It must be able to get over a 
determined value of mechanical power in heat, 
arousing the increase of the temperature of MR fluid. 
Determination of the temperature distribution is very 
important, given that MR fluid works under a limited 
temperature range. 
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FIGURE 8 BRAKING TORQUE VS NI FOR EACH ACTIVE 
SURFACE (SECOND OPTIMIZATION). 
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FIGURE 9 ELECTROMECHANICAL CHARACTERISTIC 


Considerations on the thermal equilibrium of the 
system are useful for a better clearness. 

The heat exits from the system through a convection 
exchange, while the energy enters under mechanical 
form, for effect of the sliding between the rotary parts, 
or under electrical form, for Joule effect due to the 
current in the coil. Under these conditions, the heat 
balance is: 


Q-w = — (io) 

at 

where Q is the heat quantity transferred outside in 
the time unit, W is the energy in the unit time and 
dU / dt is the variation of the inner energy of the 
system. One has: 
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Q = ~h-A,-[0(t)-0 o \ ( 11 ) 

h is the convection coefficient, 0(0 and 0o represent 
respectively the temperature in function of the time 
and the temperature of environment, while As is the 
outer surface affected by the convection flow. The 
power entering in the system can be expressed as: 

^ = ~^MR ~ 

Pmr is the mechanical power dissipated by the brake 
and Pei represents the electrical power supplied by the 
winding. The variation of inner energy of the system 
can be written as: 


dU 

dt 


d6 v 

/ m r 

dt ^ 1 


p,i 


(13) 


dd/dt is the temperature gradient, while 

^ 7Yl i * C p . represents the sum of the thermal masses 

of the system that is the product between the mass of 
every component and own specific heat. The heat 
balance of the brake is: 


-k-A,-w)-e„]HP m +p*J= 

iW „ (14) 

= > m. -c ni 

dt ^ 1 P ’ 1 

That is a differential equation of the first order at 
variable coefficients. 

1) Determination of the heat flux in the braking 

The heat generated during braking is due to the 
dissipation of both the mechanical power and the 
electrical power supplied by the winding for Joule 
effect. 

The appraisal of the heat flux generated during the 
relative motion between the surfaces of the rotor 
and of the stator in contact with MR fluid, can be 
determined calculating the variation of kinetic 
energy of the vehicle under braking. Based on the 
characteristics of the vehicle, assuming initial speed 
equal to V = 120 km/h (33,3 m/s), kinetic energy 
variation in braking until the stop is: 

(15) 

ka the inertia coefficient that takes in account also 
the inertia of the rotary mass is a function of the 
transmission ratio. Moreover the percentage of 
braking for every front brake may be calculated 
taking into account the distribution of braking force 
on the axes. The power in braking can be calculated 


depending on the time of braking. At last, the 
power transformed on heat for Joule effect in the 
coil has to be added: 

^el ~ Pcu ' Aax ' K ( 16 ) 

Vb is the volume of the coil obtained by the 
optimization, while p Cu =1,7*10 8 f 1m is the 

electric resistivity of the copper at 20 °C and Jmax is 
the current density. Contribution of this term is 
negligible in comparison with the previous. 

TABLE 5 PHYSICAL PROPERTIES OF MR140CG LORD CORPORATION. 


Base Fluid 

Oil hydrocarbon 

Color 

Dark gray 

Working temperature (°C) 

-40 - 130 

Viscosity (Pa-s) 

0,280 

Density (kg/m 3 ) 

3540 

Weight percentage of particles 

85,44% 


2) FE thermal analysis 

Thermal analysis of the brake is implemented 
using the geometry obtained from the optimization 
process. The software demands the insertion of the 
thermal properties of the materials. The properties 
of the MR fluid are reported in Table 5, while Table 
6 reports the properties of other materials. 

The loads and the boundary conditions are the 
power dissipated during the braking and the 
convective exchanges between the external surface 
of the brake and the environment. The dissipated 
power is introduced into Ansys like specific power 
in correspondence with the active surfaces of MR 
fluid. The simplificative hypothesis can be made 
that the heat generated for unity of surface is 
uniformly distributed, although the value of 
tangential stress in the fluid and the correspondent 
development of heat is variable (Cengel, 2007; 
Wiehe, Noack and Maas, 2009). In effects, the 
thermal flow is proportional to tangential stress 
that is a function of the intensity of magnetic field, 
and is proportional to the relative speed between 
the surfaces, so that the variables are numerous. A 
more accurate calculation that takes in account the 
variation of all the quantities, can be a successive 
development for a more precise temperature 
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calculation 


TABLE 6 MAIN PROPERTIES OF THE BRAKE MATERIALS. 


Properties 

CIO 

AISI304 

Al T6061-T6 

Copper 
(of the 
coil) 

Composition 

c 

0,08- 

0,13% 

Fe 

99,1- 

99,6% 

Mn 

0,3- 

0,6% 

C 0,08% 
Cr 18- 
20% 
Fe 66- 
74% 
Mn 2% 
Ni 8- 
10,5% 

A1 95,8- 
98,6% 

Cr 0,04- 
0,35% 

Cu 0,15- 
0,4% 

Fe 0,7% 
Mg 0,8-0,2% 


Density 

(kg/m 3 ) 

7860 

8000 

2700 

8920 

Yield 

strength 

(MPa) 

305 

215 

276 

48 

Ultimate 

strenght 

(MPa) 

365 

505 

310 

216 

Elastic 

modulus 

(GPa) 

205 

193 

68,9 

117 

Poisson ratio 

0.29 

0.29 

0,33 

0,355 

Thermal 

Conductivity 

(W/m-K) 

49,8 

16,2 

167 

140 

Specific heat 
(J/Kg-°C) 

448 

500 

896 

385 


Usually it is difficult to determine the convection 
coefficient, being depending on the dimensions of 
the surface of heat exchange, the speed of the air 
and its density, dynamic viscosity, thermal 
conductivity and specific heat. The existing 
relations for the determination of the value of h are 
various, based on the calculation of characteristic 
numbers in dimensionless form. The condition can 
be assumed for the brake under investigation that 
the heat exchange is of forced type and the speed of 
the air at the external surface is equal to that of 
translation of the vehicle. Convection coefficient 
can be determined by the following relationship: 


h = 


IdJ 


•Nu 


(17) 


where A is the thermal conductivity of the air equal 
to A = 0,0253 W/m°C at 20°C and Nu is the number 
of Nusselt, that is function of both Reynolds 
variable and Prandtl number Pr=0,714 (20°C). 
Density and viscosity of the air are respectively, at 
20°C: 


p a = l,224^| ?7a=L8-10- 5 ^ 
m m • s 


(18) 


Relationship (18) referred to the case of a cylinder 
having circular section, is useful for this type of 
brake; in addition Figure 10 shows the trend of 
convention coefficient versus the vehicle speed. 



FIGURE 10 CONVECTION COEFFICIENT VS VEHICLE 
SPEED. 


3) Dynamic analysis 

Typical test, in order to estimate the efficiency of 
the brakes of traditional type, is taken in 
consideration to simulate the real working of the 
brake (Aleksendric, Duboka, Gotowicki, Nigrelli, and 
Virzi Mariotti, 2006). Fifteen brake applications with 
a speed of 120 km/h till the full stop of the vehicle, 
are executed with the constant deceleration lal 
=0,75 g and a braking time tf =4.53 s and every 
single stop is followed by an acceleration period of 
the duration 28 s till a speed 120 km/h. The 
employed time in order to complete each of the 15 
cycles is equal to t cy cie=32.53 s. Table 7 shows the 
detailed list of the test for every cycle. The total 
time of the test is ttest=488 s = 8.13 min, while the 
covered distance is stest = 8129,7 m. 


TABLE 7 DETAILED LISTS OF THE BRAKING TEST FOR SINGLE CYCLE. 


Time [s] 

Speed [km/h] 

Space [m] 

0 

0 

0 

28 

120 

466,48 

28+4,53=32, 

53 

0 

541,98 


Ansys dynamic analysis demands the load 
definition and the boundary conditions in function 
of the time. It is necessary to subdivide the cargo 
history in successive loadsteps that the software 
resolves in sequential way to determine the 
temperature distribution in the brake during the 
entire test. 

During the phase of acceleration, in the hypothesis 
to neglect the friction and the generated heat due to 
the viscosity of the fluid, the brake is considered 
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subjected to convection exchanges with the 
environment. 

System is subjected to the heat load during the 
braking phase, due to the power dissipated in 
correspondence of the active surfaces and to the 
convection exchanges with the surrounding, 
having temperature 0o=2O°C. Figure 11 shows the 
course of the maximum temperature in the brake 
versus the time for the entire test. 



time [s] 


FIGURE 11 MAXIMUM TEMPERATURE IN THE BRAKE 
DURING BRAKING TEST. 

Temperature ascends quickly during the braking 
phase for every braking cycle, while declining in 
correspondence of the phase of acceleration, the 
brake disperses heat for convection. The final 
temperature to the end of the 15 cycles is equal 
0/m=125,5°C, which is lower than the value of the 
temperature limit of working of the MR fluid 
6max=150°C. Figure 12 shows the map of the 
distribution of the temperature in the brake to the 
end of the test. The maximum temperature in the 
brake at the end of the test is in correspondence to 
the inner layers of MR fluid. This is mainly due to 
the difficulty of loss of heat in such zones. Too 
much high temperatures could be obtained in the 
fluid, such as to compromise the correct operation, 
in more severe conditions of exercise (greater 
number of cycles and more elevated speeds). A 
cooling system that makes easy the removal of 
produced heat has fundamental importance. 

Conclusions 

The optimization procedure of MR brake has been 
carried out in order to verify its ability to satisfy both 
braking torque and mass requirements for a vehicle 
brake. The obtained results of thermal analysis 
highlight a maximum temperature within the limit 
operating value of the fluid. However, cooling of the 


brake to remove the heat over the active surfaces is 
necessary in the hypothesis of employment of the 
brake in more severe conditions of application. This 
could be put into effect to achieve fins in the stator, 
that increase the surface of thermal exchange with the 
environment, or a cooling of " active" type, taking 
advantage of the system already presented on the 
vehicles. 


NODAL SOLUTION 


AN SYS 


STEP=30 

SUB =5 

TIME=487 . 95 

TEMP (AVG) 

RSYS=0 

SMN =71.518 



FIGURE 12 TEMPERATURE DISTRIBUTION AT THE END OF 
THE BRAKING TEST. 


Advantages derived from the use of MR brake in place 
of a conventional brake (disc or drum) for motor 
vehicle are significant. In the first instance, MR brake 
provides reduced response time compared to 
hydraulically and/or air operated brake and it takes 
some time to raise the line pressure of the fluid to 
initiate braking. The greater response time brings 
reduction of the space, making, therefore, braking 
results more effective. The substitution of the 
conventional hydraulic braking system with an 
electrically operated system offers advantages in terms 
of simplicity of configuration and better solutions of 
housing. 

Possible future developments of MR brake deal with 
the choice of the bearings to be used between rotor 
and stator, and a complete design of the devices of seal 
of the fluid, in order to avoid potential losses, without 
compromising braking efficiency. It is moreover 
important to emphasize that the seals have not to be 
attacked by the magnetic flux, in order to avoid the 
solidification of MR fluid in their proximity with 
consequent damaging. 

Future developments of this work are: 


■ the study of the reduction of motion 
resistance due to the viscosity of MR fluid 
under normal operation, implementation 
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of a gap of variable thickness introducing 
spider legs or other. The results can be 
improved by making use of CFD. 

■ set up an analytical model or tool to 
predict thermo-dynamic behaviours of MR 
brake and build a working prototype to 
verify the performances. 
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